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Abstract 

The linked-cluster expansion technique for the high-temperature expansion of spin 
model is reviewed. A new algorithm for the computation of three-point and higher Green's 
functions is presented. Series are computed for all components of two-point Green's 
functions for a generalized 3D Ising model, to 25th order on the bcc lattice and to 23rd 
order on the sc lattice. Series for zero-momentum four-, six-, and eight-point functions 
are computed to 21st, 19th, and 17th order respectively on the bcc lattice. 

1 Introduction 

The high-temperature (strong-coupling) series expansion is one of the most successful tools for 
the study of physical systems near a critical point. 

High-temperature series are analytic; the radius of convergence is usually quite large, often 
reaching the boundary of the high-temperature phase. This property allows the application 
of powerful techniques of resummation and analytical continuation PJ, which can yield very 
precise and reliable results, provided that long series are available. It is therefore worthwhile to 
push the computation of high-temperature series as far as our algorithms and computers allow. 

The most successful technique for the computation of high-temperature series of 3D spin 
models is the linked-cluster expansion (LCE), which is well suited for the fully computerized 
approach required to reach very high orders of the expansion. 

Several detailed discussion of the LCE appeared in the literature. Wortis' review || covers 
most of the basic topics, and provides many graphical rules fit for algorithmic implementation. 
Nickel performed a remarkable computation for the 3-D Ising model on the bcc lattice ||, 
which had not been surpassed until the present work; in a more detailed paper, Nickel and 
Rehr also present several clever algorithms which we found very useful [|J . Liischer and Weisz, 
describing their application of the LCE to lattice field theory, also provide several important 
implementation hints 0. 

Unfortunately, the notation found in the literature is by no means uniform. Therefore 
we will review the relevant aspects of the LCE, which can be found in Ref. 0, not only to 
make our paper more self-contained, but also to explain notations carefully and to remark the 
correspondence with Refs. [0, |, We will follow the notations of Ref. || whenever possible. 

The paper is organized as follows: 

Sect. |2] introduces the relevant graph theory concepts and definitions. 
Sect. |3] presents the generalized Ising model we focus on. 



Sects. ||, [5], and |] review the LCE, with special focus on the two-point Green's functions. 
Sect. [F] describes our algorithm for the computation of three-point and higher Green's func- 
tions. 

Sect. |8] is devoted to programming details. 

Sect. |9] displays a (small) selection of the series generated. 

Forthcoming papers will be devoted to the analysis of the series, using the techniques pre- 
sented in Ref. , and to the generation and analysis of the series for XY systems. 

We will not give proofs of our formulae. The only nontrivial step in the proofs of Sects. |5|-[7| 
is to show that the symmetry factors compensate exactly the different number of contributions 
that may appear on the two sides of the equations; it is typically a straightforward, if tedious, 
exercise in combinatorics. The proofs of Sects. |]-|6| are given or sketched in Ref. ||. The proofs 
of Sect. [7| are especially easy, since the symmetry factor of a 1-irreducible tree graph is always 
1. 

2 Graphology 

In this section we introduce a number of graph theory concepts relevant for the LCE. We refer 
the reader to Ref. [|7| for a comprehensive introduction to the subject. 

A graph is a set of vertices and edges (also named links or bonds in the literature). Each edge 
I is incident with two distinct vertices, its extrema (we do not allow the extrema to coincide); 
the set of extrema will be denoted by dl; we write dl = {i(l), f(l)}; the choice of an "initial" 
and a "final" vertex is arbitrary. Two vertices are adjacent if they are the extrema of the same 
edge. We will denote the number of vertices and edges of a graph by v and e respectively. We 
also consider arcs or oriented edges, incident out of the initial vertex i(l) into the final vertex 
f(l). 

The valence n(i) of a vertex i is the number of edges incident with i. 

An r-rooted graph is a graph with v vertices, v > r: r roots or external vertices and v — r 
internal vertices. We will assign the indices 1, r to the roots and the indices r + 1, v to 
the internal vertices. In the drawings, roots will be denoted by open dots and internal vertices 
by filled dots. 

Two r-rooted graphs are isomorphic if there exists a one-to-one correspondence 7r of their 
internal vertices and edges such that the incidence relations are preserved, i.e., d(ir(l)) = 
{n(i(l)), 7r(/(/))} (vr(i) = i for the roots). From now on, we will identify isomorphic graph, and 
silently assume that all sets of graphs we define contain only non-isomorphic graphs. 

The symmetry factor S(Q) of an r-rooted graph Q is the number of isomorphisms of Q into 
itself, i.e. the number of of permutations of internal vertices and edges preserving the incidence 
relations. 

We will also consider p-ordered r-rooted graphs, p < r, i.e. the classes of rooted graphs 
isomorphic up to permutations of r — p roots. We will assign indices 1, p to the fixed 
roots and indices p + 1, r to the roots which can be permuted. The symmetry factor S(Q) is 
defined as the number of isomorphisms of Q into itself, including permutation of the roots p+ 1, 

r. The symmetry factor divided by (r — p)\ is called the modified symmetry factor Se{Q) 
(cfr. Ref. ||); it need not be an integer. The r-rooted graphs defined above are ordered (r- 
ordered); unless otherwise specified, we will assume that graphs are ordered. 0-ordered graphs 
are unordered. 

We will also discuss (p-ordered) r-rooted graphs whose edges and/or vertices are assigned 
a label; let us consider e.g. the case of an edge label a{l) and vertex label h{i). Two labelled 
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graphs (£?, a, b) and (£?', a', b') are equivalent (isomorphic) if there exists an isomorphism of Q 
into Q' such that a is mapped into a' and b is mapped into b' . The symmetry factor S{Q, a, b) 
is the number of isomorphisms of (Q, a, b) into itself. 

A pair of vertices % and j is connected if there exists a sequence of vertices k\, fc n , with 
k\ — i and fc n = j, and a sequence of edges Zi, / n _i such that dl a = {k a ,k a+ i}, a = 1, 
n — 1. A graph is connected if every pair of its vertices is connected. In the following, unless 
otherwise noted, we will assume that every graph is connected. 

A sequence of distinct vertices k\, k n and distinct edges h, l n is called a loop of length 
n if dl a = {k a , k a +i}, a = 1, n — 1, and <9Z n = {/c n , fci}. The number of independent loops 
(also known as cyclomatic number) of a connected graph is e — v + 1. 

A connected graph is called a tree graph if it contains no loop. A tree graph has v — 1 edges. 



3 The model 

We wish to compute the high-temperature (HT) expansion of the g-point functions of a gen- 
eralized Ising model on a D-dimensional Bravais lattice A; notice that Bravais lattices enjoy 
inversion symmetry at each lattice site. The model is defined by the generating functional 



«p(w[/i])=-n 



exp(/i i i 



ex P ( K ^2 ^3 



where 0» is a scalar field, / is an even non-negative function or distribution decreasing faster 
than exp(— 2 ) as <fi — > oo, normalized by the condition 



K = (3J, the sum runs over all pairs of nearest neighbours, and the normalization Z is fixed 
by requiring W[0) = 0. 

The connected g-point function (denoted by M. in Ref. [Q]) at zero magnetic field is defined 

by 

d q W[h] 



'11 ■■■ Ylq 



Y 



h=0 



dh h ... dhi, 



h=0 



where Xi = (xi(i), ...,xjj(i)) is the coordinate vector of the lattice site i. G q is invariant under 
the lattice symmetry group, including (discrete) translations, and under permutation of its 
arguments; it is customary to write G ? 2(xj 1 ,Xj 2 ) in the form G^^ii — Xi 2 ). We will apply the 
LCE to the computation of G q . 



4 Unrenormalized expansion 

Let us parametrize the potential / in terms of the bare vertices no(2n), defined by the generating 
function 



exp 



E 



M2n) h 2n 



d(f)f{4>) exp(/i0). 



(3) 



V (2n)! 

These quantities are named bare semi-invariants and denoted by M n ° in Ref. [0; they are named 
cumulant moments and denoted by /^n m Ref- @i A i o(2«.) = (2n — l)!!m2^, n in the notations of 
Ref. ||. Without loss of generality, we can rescale <p and K to fix /io(2) = 1. 
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For a generic r- rooted graph (?, we define the bare external vertex factor 

r 

\/ (e) (ni,...,n r ;£) = J[ Ho(n(i) + ra<), (4) 
i=i 

the bare internal vertex factor 

(£)= f[ MM*)), (5) 

j=r+l 

and the bare edge factor 

e 

L (xi, ...,x v ;Q) = Y[[K6(xi(i) -x m )], (6) 
i=i 

where 6{x) = S(\\x\\ — 1) and \\x\\ is the lattice distance between and x. 

It is convenient to focus on the contributions of r-rooted graphs to g-point functions. To 
this purpose we introduce the auxiliary r-point functions X, whose unrenormalized LCE is 

v( \ v-* vi e \n l ,...,n r ;g)V^\g)L (x l ,...,x v ;g) 
X r (x 1 ,...,x r ;n 1 ,...,n r )= ^ 2.. ^ » ( 7 ) 

where Z( r, °) is the set of all r-rooted connected graphs. 

X is invariant under simultaneous permutation of coordinates and valences: 

X r {x\, x r , n±, n r ) vY" r (x 7r (i), x^^, fi^i), n^^), 

but not over independent permutation of coordinates and valences. Furthermore, X is invariant 
over the lattice symmetry group, e.g. it is translation- invariant: 

Xi(x;n) is independent of x, and it will be denoted by X 1 (n); it will also be denoted by /i(n) 
in its role of renormalized vertex. X 2 (x±, x 2 ; rii, n 2 ) only depends on the difference x 2 — x±, 
and it will be denoted by X 2 {x 2 — Xi,ni,n 2 ). Since, by invariance under space inversion, 
X 2 (x; ni,n 2 ) = X 2 (—x; ni,n 2 ), we also have X 2 (x; ni,n 2 ) = X 2 (x; n 2 , n\). 

The sum over the location of internal vertices of the 9 functions is by definition the (free) 
lattice embedding number of Q with fixed roots: 

e 

Y[6(x i(l) - x f(l) ) = E(x u ...,x r ;g). (8) 

Therefore 

X r (xi,...,a; r ;ni,...,n r ) = 2^ ^77n ■ ( 9 ) 

Finally, the g-point functions are computed as: 

r 

Gq(xi, Xq) — Xfi^Xi^ , Xi rl , Ml, U r ) JJ^ 5 U[ {Xi n , ); (10) 

partitions i=l 



Figure 1: Example of decomposition of a graph into its 1-skeleton and 1-decorations. Open dots are 
roots, filled dots are internal vertices. 

where the g-point delta function is 



and {{in, ...,i\ Ul }, {i r ±, i rUr }} is a generic partition of {1, q} into r sets of size U\, 
u r . We will call a root bearing a factor of 5 U a u-th order root. 
The two-point function is simply 



5 Vertex-renormalized expansion 

The graph Q\i is obtained by deleting from Q the vertex i, i.e. by removing i and all edges 
incident with i. A vertex i of a rooted graph is called an articulation point if there exist 
vertices of Q\i not connected to a root. A rooted graph is called 1-irreducible if it does not 
contain any articulation point. 

Any r-rooted (r > 1) connected graph Q can be decomposed in a unique way into a 1- 
irreducible r-rooted 1-skeleton S and a 1-rooted 1- decoration for each vertex; C/ is reconstructed 
by decorating each vertex, identifying the root of its decoration with the vertex; an example is 
presented in Fig. |T[ 

Since the only 1-irreducible 1-rooted graph is the single-vertex graph, we use a different 
definition for 1-rooted graphs. A 1-rooted graph is called a 1-skeleton if it has no articulation 
points except the root. Any 1-rooted connected graph Q can be decomposed in a unique way 
into a 1-rooted 1-skeleton S and a 1-rooted 1- decoration for each vertex except the root, which 
is left undecorated. 

A 1-rooted connected graph is called a 1-insertion if Q\l is connected. 

The LCE can be reorganized by summing together all contributions from graphs having the 
same 1-skeleton, incorporating 1-decorations into renormalized vertices fi(n) = Xi(n) (named 
semi-invariants and denoted by M n in Ref. ||; fi(n) — (n — l)!!m n in the notations of Ref. 
0). The unrenormalized LCE of n(n) is given by Eq. (|7]). 

The r-point function can be computed restricting the sum in Eq. (|7]) or (|^) to the (much 
smaller) set X^ r ' 1 ' of 1-irreducible r-rooted graphs: 




1=2 



G 2 {x) = X 2 (x;l,l) + S(x)X l (2). 



X r (.?.' i . x rj 77>i, n r ^j 



(jgjO.i) x r+ i,...,x v 



E E 



Vi e) (ni, ...,n r ; Q) Vi\Q) L (xi, ...,x v ; Q) 

S(G) 
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£ n (£) K, - , n r - Q) V}*\g) K<® E(x u x r - g) 
where the internal and external renormalized vertex factors are 

r 

V} e) {n u n r - Q) = J] ^(n(i) + rw) (12) 
i=i 

and 

V?{G)= fl A*(n(i)). (13) 

i=r+l 

Eq. (0) requires a sum over all connected graphs, and therefore it is impractical for the 
computation of fi(n) at large orders of the LCE. We introduce the renormalized moments q(n) 
(named self-fields and denoted by G n in Ref. ||), defined by 

*»)= E E y °'"' g > L °';- r -^ g ) , (M) 



n .in 



where Z^ 1 -^ is the set of 1-insertions with root of valence n. The following equations hold: 

«(»)= e e ^'" (g)L t;r^° ;g) - < is » 

c e2 -(i,i) X2,...,x v ^{y) 

" n,in 

where = X^ i '°' ) flX^ 1 ' 1 ^, and X^ 1 ' 1 ^ is the set of 1-rooted 1-skeletons; 

OO 1 oo oo 

H(n) = fio(n) + — J2 ■■■ J2 ■■■ Voin + h + ■■■ + l a ). (16) 

S=l S - h=l l s =l 

Since q{2n — 1) = 0, and odd values of n and k do not contribute to Eq. fllEp. g and /x can now 
be computed recursively in parallel order by order in K, since, once Eq. (^) is expanded in 
powers of K and truncated, the coefficient of the highest power of K of q in the l.h.s. depends 
only on coefficients of lower powers of K of \x in the r.h.s. 



6 Edge-renormalized expansion 

A pair of distinct vertices % and j of a rooted graph Q is called an articulation pair if there 
exist vertices of Q\i\j not connected to a root, or if i and j are joined by more than one edge. 
A rooted graph is called 2-irreducible if it does not contain any articulation pair. 

Any 1-irreducible r-rooted (r > 2) graph Q can be decomposed in a unique way into a 
2-irreducible r-rooted 2- skeleton S and a 1-irreducible 2-rooted 2- decoration for each edge (ori- 
ented in a canonical way, e.g. by choosing i(l) < /(/)); Q is reconstructed by replacing each 
edge with its decoration, identifying the first and second decoration root with the initial and 
final vertex of the edge respectively. An example is shown in Fig. 0. 

We use a different definition for 2-rooted graphs, since the only 2-irreducible 2-rooted graph 
is the bond graph (no internal vertices and only one edge); the roots of all other 1-irreducible 
2-rooted graphs are an articulation pair. We call a 1-irreducible 2-rooted graph a 2-skeleton 
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Figure 2: Example of decomposition of a graph into its 2-skeleton and 2-decorations. Open dots are 
roots, filled dots are internal vertices. 

if it does not contain any articulation pair except the pair consisting of the two roots. Any 
1-irreducible 2-rooted graph can be decomposed in a unique way into a 2-rooted 2-skeleton S 
and a 1-irreducible 2-rooted 2- decoration for each edge except edges connecting the roots, which 
are left undecorated. 

The LCE can be reorganized by summing together all contributions from graphs having the 
same 2-skeleton, incorporating all 2-decorations into renormalized edges. 
We start by decomposing Eq. ( |PTD for r > 1 into 

w/ / n Vi\Q) L Q {x 1 ,...,x v ;g) 
W r {x 1 ,...,x r ;n 1 ,...,n r ) = E FTTJn » ( 17 ) 

c T (r,l) X r+l ,...,X v ^{y) 

y ^-^n^ , . . . ,nr 



Sl,...,S r 



YlfJ,(rii + Si 



.i=i 



W r (x u ...,x r ;si, ...,s r ), (18) 



where T%f) Ur is the set of 1-irreducible r-rooted graphs with roots of valence ni, n r . W r 
enjoys the same symmetry properties of X r . Wi(x,n) is undefined. 

W r can be computed by assigning an initial and a final valence ii, fi to each oriented edge 
of a 2-rooted graph; the valence i\ is incident with i{l) and fi is incident with f(l): 

r 

W r (x 1 ,...,x r ;n 1 ,...,n r )= ^ ^ JJ 5(n; - Ui(i u i e , / x , / B ; £)) 

gex( r ' 2 ) ZV+lv^uSH,— ,ie,/l,— ,/e *=1 

Vj^gi, Ze, A, .-, / e ; g) L 2 (xi, x v ; i u ...,i e , f u f e ; g) 

S(f?) ' ( j 

where X^' 2 ^ is the set of 2-irreducible r-rooted graphs, Ui is the sum of all the valences incident 
with the vertex i, 

V 

Vj % \ii,...,i e ,fi,...,fe]G)= II KMhi-^eJu-JejS)), (20) 

i=r+l 

and 

e 

L 2 (xi, ...,x v ;i u ...,i e , f u f e ;Q) = Y[W 2 (x ii i ) -x f (iy,ii,fi). (21) 

i=i 
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For r = 2 Eqs. (|19D and (21) require a slight modification: we define Z^ 2 ' 2 ) as the set of 
2-rooted 2-skeletons; for edges incident with at least one root, we replace W 2 {x\ — x 2 ; h f) with 
K6(xi - x 2 )5n5f 1 . 

The sum over graphs can be restricted to a subset of Z( r ' 2 ); we will discuss here the case 
r = 2; the next Section will be devoted to the case r > 3. Let us start by classifying 1-irreducible 
2-rooted graphs into several classes. 

An internal vertex i of a 1-irreducible 2-rooted graph Q is called a nodal point if the roots 
of Q\i are not connected. A 1-irreducible 2-rooted graph is nodal (also named articulated or 
separable in the literature) if it contains one or more nodal points; otherwise it is non-nodal. 

A 1-irreducible 2-rooted graph is simple if {?\1\2 is connected, and 1 is not adjacent to 2. 
By definition, all nodal graphs are simple. A 1-irreducible 2-rooted graph is a ladder graph if 
it is not simple and it is not the bond graph. 

A 1-irreducible 2-rooted graph is elementary if it is both simple and non-nodal. 

We have divided 1-irreducible 2-rooted graphs into four disjoint classes: bond, nodal, ladder, 
and elementary graphs. Let us separate the contributions to W 2 according to the four classes: 

W 2 (x; ni ,n 2 ) = W*°{x; m, n 2 ) + W 2 no (x; n x , n 2 ) + W^{x; n x ,n 2 ) + Wf(x; n u n 2 ), 

i.e. bond, nodal, ladder, and elementary contributions respectively. 

The bond contribution is trivial. Nodal contributions can be factorized into a product of 
non-nodal contributions: 

W™(x;n u n 2 ) = ]T W^(x 3 ;m, k) + k) W™(x - x 3 ; i 2 , n 2 ) 
+ E W^{x 3 -n x ,k)^{k + i 2 ) 
x W* n (x A - x 3 ; i 2 , z 3 ) M*3 + k) W™(x - x A - z 4j n 2 ) + ... , (22) 
which can be written recursively as 

W^°(x;n 1 ,n 2 ) = ]T ^ 2 nn (x 3 ; m, k) ii{k + i 2 ) W 2 (x - x 3 ;* 2 , n 2 ). (23) 

Likewise, ladder contributions can be factorized into a product of non-ladder contributions: 

00 1 s 

Wf(x; ni ,n 2 ) = J2- E S( ni -U=ih)6(n 2 -j:Lift) II ^0*. ( 24 ) 

s=2 S - h,...,i s ,h r ..,fs t=X 

Elementary contributions can be computed by setting r = 2 into Eq. (|1J), and restricting 

(2 2) 

the sum to ' , the set of elementary 2-rooted 2-skeletons. The sum can be further restricted 
to , the set of unordered elementary 2-rooted 2-skeletons, provided that we replace S(Q) 
with Se{Q) and we symmetrize the result: 

1 2 
Wf(x;n u n 2 ) = - E E ]j5(n t - Ui(i u ...,i e J u ...J e ;G)) 

a —(2,2) x 3 ,...,X v ;i 1 ,...,i e ,fl,...,fe '=1 
yt 0,el 

V 2 (l \k, ie, fl, ■■; fe] Q) L 2 (x h X v ] k, ■■; le, fl, ■■; fe] Q) . , OK \ 

x ^ +n 1 -n, (25) 
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Figure 3: Decomposition of a 2-irreducible graph into a compact kernel and 2-irreducible tree graphs. 



The last ingredient we need for a fully edge-renormalized expansion are the renormalized 
vertices /x(n); they can be computed by combining Eq. (|16|) with 



q(n 1 +n 2 ) 



n\ \ n 2 ! 



r (WT(0;ni,n 2 ) + W 2 el (0;n 1 ,n 2 )) 



(26) 



(ni + n 2 ) 

reflecting the fact that the contributions to q{ni+n 2 ) in Eq. ( |Hf) can be obtained from the 
contributions to W^O; 711,712) in Eq. (|T/]) by identifying the roots of the 2-rooted graph, i.e. 
by suppressing the second root and reattaching all the edges incident with it to the first root, 
provided that the roots are not adjacent. 

Expanding in powers of K Eqs. (|23|) , (|2^) , (]25|) , (|26|) , and (0), we can compute Wf^x; n 1; n 2 ), 



W / 2 la (x;ni,n 2 ), <?(n), 



which involves a summation over graphs is Eq. (25), where the sum only runs over Xq 
relatively small set. 



and /i(n) in parallel order by order in K. The only step 

(2,2) 
el 



a 



7 Three-point and higher functions 

A vertex (internal or external) i of a 1-irreducible r-rooted graph Q is called a nodal point if 
is not connected. By definition of 1-irreducibility, each connected component of Q\i must 
contain at least one root. A 1-irreducible r-rooted graph Q is nodal if it contains one or more 
nodal points; otherwise it is non-nodal. 

In the rest of this section, we will assume that every graph is 2-irreducible. 

A nodal point j of a 2-irreducible r-rooted graph Q is called a tree-insertion point if at least 
one of the connected components of Q\j is a tree graph. The order t of j is the number of roots 
of Q contained in all tree graph components of plus 1 if j itself is a root. A 2-irreducible 
r-rooted graph is called compact if it contains no tree-insertion points. 

Let us consider a 2-irreducible r-rooted graph Q . We generate a compact 2-irreducible 
graph Q f , called the compact kernel of Q, by removing all the tree graphs attached to every 
tree- insertion point, and promoting all internal tree-insertion points to root. Q is obtained by 
attaching a 2-irreducible tree graph to each root of Q' . If Q' is not a tree graph (and therefore 
it has at least 3 roots), the decomposition is unique; otherwise, Q itself is a tree graph. An 
example is shown in Fig. 0. 

By summing all contribution of 2-irreducible graphs with the same compact kernel, we can 
write W q as 

W q (xi, ...,x q ; m, n q ) = Wq(x!, x q ] n 1; ...,n q ) 

r 

+ Yl W^ (y 1 ,...,y r ;i 1 ,...,i r )]J Y Ul+1 (yi,x in , ...,x ilui ; 11,11^, ...,n hui ), 

partitions yi,...,y r ,ii,...,i r 1=1 
r>2 
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Figure 4: Graphical representation of Eq. (29) for q = 5. 



(27) 

where W q T (xi, x q ; ni, n q ) is the tree graph contribution to W q , W^°(yx, y r ; i\, i r ) is 
the compact graph contribution to W r , and 

Y t+1 (y,x h ...,x t ;i,nx, ...,n t ) = ^^{i + j) W£ a (y,zi, ...,x t ; j,m, ...,n t ) 

j 

+ Y. 6 (y - x 6 ( n ± -i-j) w t T ( x u W,n 2 , ...,nt), 
j 

Y 2 (y,x;i,n) = £ fi(i + j) W 2 tT (y, x; j,n) + 5(y - x) 5(n - i) (28) 

3 

(the first and second term correspond to an internal and external i-tree-insertion point respec- 
tively); notice that W 2 = W 2 . Eq. fl2"T| ) can be combined with Eq. (|II]) to give 

Gq{X\ 1 Xq) (jr (Xl, Xq) 

r 

+ E ^ (y 1} ...,y r ;ii,...,i r )Y[Z Ul+1 (y h x ill ,... } x ilui ;ii) } (29) 

partitions yi,...,y r ,ii,---,i r 1=1 
r>2 

where 

Z t+ x(y,x 1 ,...,x t ;i) = E E ^H-iG/, x in , x irl ; i, s u s r ) 

partitions si,...,s r 
of {l,...,i} 
r 

x n^^+^^Ki)-^^)- ( 3 °) 

Z t+1 (y,Xi, ...,x t ;i) is symmetric under permutations of Xi,...,a; t and lattice symmetries, e.g. 
simultaneous translation of y and Xx, x t . 

These formulae can be written graphically, according to the rules presented in Table |I| A 
sum over all dummy coordinates y and all dummy valences i, f and a sum over all inequivalent 
permutation of external coordinates x or coordinate-valence pairs x, n are understood. Notice 
that, despite the graphical notation, all pairs of roots of W q n and W q ° are equivalent. 

Eq. ( p9|) can be expressed by writing all polygons with a letter "c" having 3 to q vertices, 
placing a crossed dot with a positive integer label at each vertex in all inequivalent ways, the 
sum of the labels being q, and adding the tree contribution. The case q = 5 is shown in Fig. |J 

G q and Z q can be computed by adding the contributions of all 2-irreducible r-rooted tree 
graphs with roots labelled by positive integers with sum g; for Z q , the first root must be drawn 
as a square. The case Z 3 is shown in Fig. |. 
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Table 1: Graphical rules, v is the sum of all the valences incident with the vertex. The argument Z{ 
of Wq n is the coordinate x or y associated with the symbol placed at the vertex i of the polygon; the 
argument y^ of WI 20 is the variable y of the function Z q+ \ placed at the vertex i. 
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comment 


variables 


factor 
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internal vertex 


y 











unlabelled root 


x, n 


S(n — u) 




O 
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Figure 6: Compact contributions to W£° and corresponding nodal skeletons. 



The next step is to write an expression of Wq° in terms of W™. For q = 3 we have simply 
= W™. The case q = 4 is shown in Fig. |6|. For larger values of q, the number of non-nodal 
contributions to grows rapidly, and a systematic approach is needed. 

Let us define for a (connected or non-connected) graph Q and a vertex i the graph Q/i: let 
Q* be the connected component of Q containing i\ if Q*\i is connected, set Q/i = Q; otherwise, 
for each connected component Q% of generate the graph Q k by adding a new vertex, internal 
or external like i, and joining it to all the vertices adjacent to i in Q*, by the same number 
of edges; replace Q* with the connected components Qk- The edges and vertices of Q/i are in 
one-to-one correspondence with the edges and vertices of i, except for the new vertices which 
all correspond to i (a nodal point of Q). Notice that Q/i/j = Q/j/i- 

Let Q be a compact 2-irreducible r-rooted graph (r > 3) with t nodal points ii,...,i t . 
Observe that all the connected components Q\ of Q = Qji\j...ji t are non-nodal. Generate a 
new graph T, the nodal skeleton of Q, in the following way: for each Q\ with v > 3, containing 
n roots corresponding to non-nodal roots of Q, write a root / of T with a label n > 0; for 
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Figure 7: A compact graph, its non-nodal components, its nodal skeleton, and the contribution to 
W co of the set of graphs sharing the nodal skeleton. 




2122 122 1 2 



Figure 8: Nodal skeletons contributing to W§°. 

each node ik of Q write an unlabelled vertex ik of T, internal or external like i^; join ik to all 
the labelled roots / such that Q\ contain a vertex corresponding to if., and with all unlabelled 
vertices which are adjacent to ik in Q. 

A nodal skeleton is a connected 1-irreducible tree graph, but it is not in general 2-irreducible 
(it may contain 2-valent internal vertices). A nodal skeleton enjoys the following properties: 
each 2-valent internal vertex is adjacent to a labelled root; labelled roots are never adjacent; 
unlabelled roots are at least 2-valent; m-valent roots with label n satisfy n + m > 3. Every 
nodal skeleton can be generated by adding labels to some of the roots and by splicing 2-valent 
internal vertices into a 2-irreducible tree graph. 

Every 1-irreducible tree graph, with some roots carrying a non-negative integer label, sat- 
isfying the above properties, is the nodal skeleton of a set of g-rooted 2-irreducible compact 
graphs, with q equal to the sum of the labels plus the number of unlabelled roots. The con- 
tribution of this set to can be computed by replacing each m-valent root of T carrying a 
label n with an (n + m)-sided polygon whose vertices are the m vertices adjacent to the root 
and n new (unlabelled) roots, and applying the rules of Table [I]. 

An example of the construction of the nodal skeleton and its evaluation is presented in Fig. 
[7|. The set of all nodal skeletons contributing to W£° is shown in Fig. [| see also Fig. || 
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We could carry further the reduction of the set of graphs to be summed over, e.g. by 
identifying ladder graphs along the lines of Sect. This is rather complicated for arbitrary 
r, and goes beyond the scope of the present work. Moreover, the zero-momentum projection 
described below is not applicable to the ladder graph reduction. Therefore we compute 
by restricting the sum of Eq. ([19]) to 2 \ the set of non-nodal 2- irreducible r-rooted graphs. 

The above considerations can be simplified considerably if we are only interested in moments 
of the g-point functions, e.g. 

Xq = ^gO^lj •••) x q)i 

M 2 (9) = (x 1 -x 2 ) 2 G q (x 1 ,...,x g ). (31) 

X2,—,Xq 

Let us introduce the moments of Z q , W q ° and W q n : 

C q ( n ) = E Z q (x u ... 1 x q ;n), 
<K,-,n,) = E W q M (x 1 ,...,x q ;n 1 ,...,n q ), (32) 

X2,—,X q 

the corresponding definition for x\ r an d wj" 1 ^!, ■■■,n q ), all independent second moments, etc. 
Eq. fl29D can be projected over zero momentum to give 

X fl = X? + E E T^—( E ^°(i 1: ...,i r )f[Cu l+ i(ii). (33) 

r =3 «!<...<«,. 11/=1 u l- u,...,i r 1=1 

ui + ... + u r — q 

The computation of ( q , x q r , an d ui q ° is also easy; the graphical rules can be immediately 
projected over zero momentum, suppressing all coordinates x and y and removing all space delta 
functions; the only nontrivial part is the counting of the number of inequivalent permutations 
of roots. 

C q (n) can be computed by summing over all inequivalent 2-irreducible 1-ordered r-rooted 
tree graphs, with the roots labelled by positive integers u\, ...,u r with sum q. The number of 
inequivalent permutations of the roots 2, r is 

where S is the symmetry factor of the labelled graph. 

The computation of x q T { n ) is very similar, but we sum over unordered graphs, and the 
number of inequivalent permutations of the roots is 

gj 

u) q °(ni, ...,n q ) can be computed by summing over all inequivalent unordered r-rooted nodal 
skeletons, with p roots labelled by non- negative integers U\, u p with r — p + u\ + ... + u p = q, 
with a weight 1/S, where S is the symmetry factor of the labelled graph, with unlabelled roots 
assigned an arbitrary distinct label (e.g. —1). 

Finally, we can compute u^ n (ni, ...,n q ) by summing over unordered non-nodal 2-irreducible 
g-rooted graphs, provided that we use the modified symmetry factor and we symmetrize the 
result under permutation of the valences. 
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The computation of the second moment of the above quantities proceeds along the same 
lines. The factor (x{ — Xj) 2 in Eq. fl3~T| ) is dealt with in the following way: the two roots are 
connected by a chain of terms with a space structure of the form 

53 ( X i ~ Xj) 2 fl(Xi - y X ) f 2 {Vl - y 2 ) ••• fn+liVn - Xj) 
Xj,yi,...,y n 

(we have dropped the dependency on coordinates lying outside the branch connecting i with 
j). Let us write 

(xi - xj) 2 = (xi - yx) 2 + (yi - y 2 ) 2 + ... + (y n - Xj) 2 + cross terms. 

The cross terms do not contribute to the sum, and the result is 

, (2) ,(o) AO) AO) A2) AO) AO) AO) A2) 

Jl J2 •■■in+1 T il J2 • • ' in+1 ' •• ■ ' J 1 J2 •••Jn+1) 

where 

yf^E /<(*), f! 2) = J2y 2 m- 

y y 

Therefore we can compute the second moment by taking each contribution to the zero-momentum 
quantity, promoting one of the zero-momentum factors along the branch connecting the roots 
% and j to second moment, and summing over all possible choices. 

By dealing with moments, we avoid the need of storing all the values of Z, W co , and 
W nn , which can rapidly exhaust all available memory. The extension to higher moments is 
straightforward but cumbersome. 

8 Programming details 

We wrote a set of computer programs to implement the automatic evaluation of the edge- 
renormalized LCE on the simple cubic lattice (sc) and on the body-centered cubic lattice (bcc) 
(3D); the same programs evaluate the LCE on two different representations of the square lattice 
(2D) and on the ID lattice. 

The computation of g-point functions is performed for a generic potential, keeping n (2n) 
symbolic; each term of the series is a polynomial in fi (2n) with rational coefficients. We also 
implemented the same computation for a specific potential; this requires much less memory and 
is somewhat faster (up to 30%), but not enough to give up the flexibility of a generic potential. 

To speed up search and insertion into ordered sets of data, graph sets and polynomials in 
/i (2n) are implemented as AVL trees (height-balanced binary trees) (cfr. e.g. Ref. 0, Chapt. 
6.2.3), using the ubiqx library. Rational numbers and (potentially) large integers are handled 
by the GNU multiprecision (gmp) library. 

Given the complexity of the procedure, it is crucial to perform a number of checks in order 
to flush out all algorithm and program errors. In ID our series are compared with exact results 
for the spin-1/2 [|9] and the spin-1 |H| Ising model; this is already a very stringent check, 
especially of the graph sets (cfr. Ref. Q). In 2D, our results are compared with the series for \ 
and M 2 for spin-1/2 published in Ref. In 3D, our results are compared with the lower-order 
series already available, for \ an d M 2 for specific potentials in Refs. || U |TlJ, and for \i M 2 , 
and Xq f° r a generic potential in Ref. |J. q(n) can be computed from different combinations of 
rii and n 2 in Eq. (p6|); their agreement is non-trivial. Finally, for the spin-1/2 Ising model on 
any lattice, the series for G q , rewritten in terms of v = tanhK, must have integer coefficients. 
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Figure 9: The simplest elementary 2-rooted 2-skeleton. 



8.1 Graph generation 

A program generates the table of all unordered elementary 2-rooted 2-skeletons contributing to 
the desired order; the algorithm follows Ref. H. Starting from the graph drawn in Fig. [| we 
apply recursively the following modified Heap rules ||12|| : 



(a) join any two distinct vertices by a new edge, provided the two vertices are not already 
adjacent; 

(b) insert a new internal vertex on any edge and join it to any vertex, excluding the edge 
extrema, by a new edge; 

(c) insert two new internal vertices on any two distinct edges, and join them by a new edge; 

(d) do not join the roots by a new edge. 

The modifications to the original Heap rules (a), (b), and (c), marked in italics, prevent the 
generation of 2-reducible graphs. Rule (d) prevents the generation of ladder graphs. 

The reduction of graphs to canonical form is performed using a generalization of the algo- 
rithm of Ref. H| . Graphs are stored in a compact form similar to the one of Ref. || . 

To reduce the proliferation of graphs at higher orders, it is extremely important to know the 
order ("strict bound") o s at which a given graph will enter in the expansion (it is not trivially e, 
since we require even valence of all internal vertices, and, being interested in bipartite lattices, 
even length of all loops). 

We must also keep in mind that some graphs don't contribute at the desired order, but 
graphs generated from them might contribute. We define the "Heap bound" oh{G) as the 
minimum of o s on the set of graphs including Q and all graphs generated from it. We also 
define the two bounds for (G,t]), i.e. the minimal order when the vertex i is forced to be 
embedded in a lattice site of parity nil). 

We apply the modified Heap rules to the graph Q in the following way: assign a parity label 
77(2) to each vertex, in all the ways compatible with the Heap bound; apply the modified Heap 
rules assigning all possible parity labels to the new vertices; discard immediately the generated 
graph-parity pairs not satisfying the Heap bound; discard vertex parity information and store 
the generated graphs not isomorph to previously generated graphs. Finally, save into a file only 
the graphs satisfying the strict bound. 

The generation of the elementary 2-rooted 2-skeletons contributing to the 25th order re- 
quired 41 hours of computation on one CPU of a Compaq ES-40, and ca. 300 Mbytes of RAM. 
The number of inequivalent elementary 2-skeletons for each order of o s is reported in Table |2|. 

A similar program generates all unordered non- nodal 2-irreducible r-rooted graphs for r > 3. 
The table is initialized by applying rules (a) and (b) recursively, starting from each 2-irreducible 
r-rooted tree graph, until the result is non- nodal. Rules (a), (b), and (c) are then applied 
recursively. The number of inequivalent non-nodal 2-irreducible r-rooted graphs for each order 
of o s is reported in Table |[ 

We remark that this graph tables can be used for the LCE of any spin model with <fi — > — <p 
symmetry on a bipartite lattice in any dimension. 

The generation of all the required families of tree graphs is straightforward. 
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Table 2: Number of inequivalent unordered elementary 2-rooted 2-skeletons (r = 2), or unordered 
non- nodal 2-irreducible r-rooted graphs (r > 2), for each order of the strict bound o s . The number of 
graphs satisfying the Heap bound is typically 20% to 30% higher. No graph in any of these sets has 
a bound lower than 4, 
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8.2 Computation of the g-point functions 

A separate program reads the table of elementary 2-rooted 2-skeletons and computes all com- 
ponents of Wi. 

The evaluation of \x and of bond, nodal, and ladder contribution to W2 is a straightforward 
application of the formulae of Sect. |6|. 

The evaluation of elementary contribution dominates the computation time, and must be 
optimized as much as possible. Assume that all lower-order contributions to W2 have been 
computed. For each unordered elementary 2-skeleton Q with o s not larger than the desired order, 
all inequivalent assignations Q, (n, I) of edge valence parity n and length parity I compatible 
with the desired order, with even length of all loops, and with even valence of all internal 
vertices, are generated. We have implemented two different algorithms for the computation of 
the contribution of Q, (n, I) to the two-point function. 

The first algorithm is essentially the one used by Nickel and Rehr in Ref . [§J : all inequivalent 
1-irreducible 2-rooted graphs with a 2-skeleton compatible with Q, (n, I) are generated, and their 
contributions are computed according to Eq. (]TT|) . In the second algorithm, the contribution is 
computed according to Eq. (|25|), and 1-irreducible graphs are not needed. The first algorithm 
is more efficient for 2-skeletons with large o s , while the second algorithm is more efficient for 
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small o s ; for each value of o s we select the algorithm which is (presumably) more efficient. On 
the sc lattice, the speed-up obtained over the use of either algorithm for all skeletons grows with 
the order, and is about a factor of 4 at order 23. On the bcc lattice, the first algorithm is very 
efficient, since the embedding number factorizes into a product of 1-dimensional embedding 
numbers ||]; we still use the second algorithm for the simplest 2-skeletons (o s < 10), since the 
computation of the corresponding 1-irreducible 2-rooted graphs contributing to orders higher 
than 21 is extremely time- and memory-consuming. 

Keeping in RAM all components of W2 for a generic potential would be problematic. Most 
of these components are needed only to compute nodal and ladder contributions, and can be 
kept on disk; keeping in RAM just the components needed to compute elementary contribution 
is manageable. 

The computation of the 25th-order LCE for the two-point function on the bcc lattice re- 
quired ca. 400 hours of computation, and ca. 700 Mbytes of RAM. 

A similar program reads the table of non-nodal 2-irreducible r-rooted graphs and the compo- 
nents of W2, and computes u™. The computation of Xq is then straightforward. We computed 
Xi, X61 an d Xs f° 21st, 19th, and 17th order respectively on the bcc lattice. 

The computation of the same quantities on the sc lattice is much slower (but does not 
requires more RAM); so far, we obtained W2 to 23th order, with an effort not much smaller 
than the 25th order on the bcc lattice. The computation of W% and Xq to the same orders as 
on the bcc lattice is in progress, but it will require a non-trivial amount of time. 



9 Selected results 

All high-temperature series computed in the present work are available for the most general 
potential, in the form of polynomials in the bare vertices no(2n). The general results are 
extremely lengthy, and are only useful for further computer processing. 

We present here a selection of high-temperature series for the spin- 1/2 Ising model, i.e. for 

/(0) = |(5(0+1) +5(0-1)). 

For sake of compactness, all the series are written in terms of v — tanhfT. Series for other 
specific potentials are available upon request from the author. 

Although we computed all components of G2(x,y), we report here only x = X2 and M2 = 
(cfr. Eq. (0)). For the g-point functions, we only computed Xq- 

On the bcc lattice, we obtained 

x = 1 + 8v + 56 v 2 + 392 v 3 + 2648 v 4 + 17864 + 118760 + 789032 

+ 5201048 v 8 + 34268104 v 9 + 224679864 v 10 + 1472595144 v 11 + 9619740648 v 12 
+ 62823141192 v 13 + 409297617672 v 14 + 2665987056200 v 15 + 17333875251192 v w 
+ 112680746646856 v 17 + 731466943653464 v 18 + 4747546469665832 v 19 
+ 30779106675700312 v 20 + 199518218638233896 v 21 + 1292141318087690824 v 22 
+ 8367300424426139624 v 23 + 54141252229349325768 v 24 

+ 350288350314921653160 v 25 + 0(w 26 ); (34) 
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M 2 = 8 v + 128 v 2 + 1416 v 3 + 13568 v 4 + 119240 v 5 + 992768 v 6 + 7948840 v 7 
+ 61865216 + 470875848 v 9 + 3521954816 v 10 + 25965652936 v 11 
+ 189180221184 v 12 + 1364489291848 v 13 + 9757802417152 v 14 
+ 69262083278152 v 15 + 488463065172736 v 16 + 3425131086090312 v 17 
+ 23896020585393152 v 18 + 165958239005454632 v 19 + 1 147904794262960384 v 20 
+ 7910579661767454248 v 21 + 54332551216709931904 v 22 
+ 372033905161237212392 v 23 + 2540342425838560175616 v 24 
+ 173014572071 10720278440 v 25 + 0(v 26 ); (35) 



-X4 = 2 + 64v + 1168t; 2 + 16576v 3 + 201232 v 4 + 2204608v 5 + 22411504 v 6 
+ 215447872 v 7 + 1981980688 v 8 + 17602809920 v 9 + 151865668752 v 10 
+ 1278888344256 V 11 + 10550227820400 v 12 + 85510907958720 v 13 
+ 682500568307184 v 14 + 5374496030148928 v 15 + 41821018545214608 v 16 
+ 321992795063663936 v 17 + 2455641803116052752 v 18 + 18567879503614668736 v 19 
+ 139310655514229882000v 20 + 1037854026688655887552 v 21 + 0(v 22 ); (36) 



Xe = 16 + 1088 v + 36416 v 2 + 853952 v 3 + 15974528 v 4 + 255491264 v 5 + 3638767040 v 6 
+ 47395195712 v 7 + 574950589568 v 8 + 6581949043264 v 9 + 71803170318144 v 10 
+ 752047497945024 V 11 + 7606707093034368 v 12 + 74649010982738112 v 13 
+ 713458387977120192 v 14 + 6661638582474716480 v 15 + 60923519621981242752 v 16 
+ 546923327751320201536 v 17 + 4828463182433394315584 v 18 

+ 41987611565592990702272v 19 + 0(v 20 ); (37) 



-Xs = 272 + 31744 v + 1673728 v 2 + 58110976 v 3 + 1538207872 v 4 + 33584739328 v 5 
+ 634387677184v 6 + 10699575811072 v 7 + 164723097021568 v 8 
+ 2352360935459840 v 9 + 31540880634427392 v 10 + 400802365468148736 v 11 
+ 4862781935250449280 v 12 + 56665753776838026240 v 13 
+ 637305912177206767104 v 14 + 6945658883867865975808 v 15 
+ 73600395257678784586368v 16 + 760476823195422275111936 v 17 + 0(v 18 ). 

(38) 

On the sc lattice, we obtained 
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x = 1 + 6 r; + 30 + 150 + 726 + 3510 1; 5 + 16710 1> 6 + 79494 v 7 + 375174 v 8 
+ 1769686 v 9 + 8306862 v 10 + 38975286 v 11 + 182265822 v 12 + 852063558 v 13 
+ 3973784886 v 14 + 18527532310 v 15 + 86228667894 v 16 + 401225368086 v 17 
+ 1864308847838 v 18 + 8660961643254 v 19 + 40190947325670 V 20 
+ 186475398518726 v 21 + 864404776466406 v 22 + 4006394107568934 v 23 
+ 0(v 24 ); (39) 



M 2 = 6 v + 72 v 2 + 582 v 3 + 4032 v 4 + 25542 v 5 + 153000 v 6 + 880422 v 7 + 4920576 v 8 

+ 26879670 V 9 + 144230088 v 10 + 762587910 v 11 + 3983525952 v 12 + 20595680694 v 13 

+ 105558845736 v 14 + 536926539990 v 15 + 2713148048256 v w + 13630071574614 v 17 

+ 68121779384520 v 18 + 338895833104998 v 19 + 1678998083744448 v 20 

+ 8287136476787862 v 21 + 40764741656730408 v 22 + 199901334823355526 v 23 

+ 0(v 24 ). (40) 
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